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ABSTRACT 

We present a general approach to the modelling of the brightness and polariza- 
tion structures of adiabatic , decelerating relativistic jets, based on the formalism of 
iMatthews fc Scheuerl |l990). We compare the predictions of adiabatic jet models with 
deep, high-resolution observations of the radio jets in the FRI radio galaxy 3C31. 
Adiabatic models require coupling between the variations of velocity, magnetic field 
and particle density. They are therefore more t ightly constrained than the models pre- 
viously presented for 3C 31 bv lLaing fc Bridle! l|2002a|) . We show that adiabatic models 
provide a poorer description of the data in two crucial respects: they cannot repro- 
duce the observed magnetic-field structures in detail, and they also predict too steep 
a brightness decline along the jets for plausible variations of the jet velocity. We find 
that the innermost regions of the jets show the strongest evidence for non-adiabatic 
behaviour, and that the adiabatic models provide progressively better descriptions of 
the jet emission at larger distances from the galactic nucleus. We briefly discuss physi- 
cal processes which might contribute to this non-adiabatic behaviour. In particular, we 
develop a parameterized description of distributed particle injection, which we fit to 
the observed total intensities. We show that particles are preferentially injected where 
bright X-ray emission is observed, and where we infer that the jets are over-pressured. 

Key words: galaxies: jets - radio continuum: galaxies - magnetic fields - polarization 
- MHD - acceleration of particles 



1 INTRODUCTION 

This paper is one of a series whose aim is to de- 
velop quantitative mod els of jets in low-luminosity (FRI; 
iFanaroff fe Rilevl 11974 radio galaxies on the assumption 
that they are decelerating relativistic flows. Our hypothesis 
is that the jets are close enough to being intrinsically sym- 
metrical, axisymmetric and antiparallel that the observed 
differences between them are dominated by the effects of rel- 
ativisti c aberration. This hypothesis is motivated by the re- 
sults of lLaing et al.l il999Tl , who presented a statistical study 
of jets in the B2 sample of radio galaxies and showed that 
the observed correlations between fractional core flux den- 
sity and side-to-side asymmetries in intensity and width are 
consistent with jet deceleration and the presence of trans- 
verse velocity gradients. They inferred that the jets slow 
from ~0.9c where they start to expand rapidly to ~0.1c over 
distances ~1 to lOkpc and that the deceleration scale is an 
increasing function of jet power. In lLaing fc Bridlel {2002a, 
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hereafter LB), we demonstrated that an intrinsically sym- 
metrical, relativistic jet model provides an excellent descrip- 
tion of the total intensity and linear polarization observed 
at 8.4 GHz from the jets of the FRI radio galaxy 3C 31. We 
were able to estimate the angle to the line of sight and the 
three-dimensional distributions of velocity, emissivity and 
magnetic-field structure. By combining this kinematic model 
with a description of the external g as density and pressure 
derived from Chandra observations iHardcastle et af]l2002l) 
and using conservation of particles, energy and momentum, 
we demonstrated that the jet deceleration could be pro- 
duced by entrainment of thermal matter and we derived the 
spati al variations of pr essure, density and entrainment rate 
l|Laing fc Bridlell2002bl) . 

The free models developed by LB for 3C 31 were de- 
signed to fit the observed images without embodying specific 
preconceptions about the (poorly known) internal physics. 
In those models, we adopted simple and arbitrary functional 
forms for the spatial variations of velocity, synchrotron emis- 
sivity and field ordering and allowed the emissivity and field 
ordering to vary separately as smooth functions of position. 
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To proceed beyond such purely empirical descriptions of the 
jets, we must make further assumptions about processes that 
affect the components of the modelled emissivity - the rel- 
ativistic particle energy spectrum and the strength of the 
magnetic field. The separation of the emissivity into its com- 
ponents is ill-determined unless inverse-Compton emission 
can be detected from the synchrotron-emitting electrons, in 
which case the particle density and field strength can be 
determined independently. Unfortunately, inverse-Compton 
emiss ion from the inne r jets in 3C 31 is too weak to detect 
jHardcastle et alJl20o3) . There are, as yet, no prescriptive 
theories for dissipative processes such as particle acceler- 
ation and field amplification or reconnection in conditions 
appropriate to FR I radio jets. 

The energy-loss processes for the radiating particles can 
be quantified, however. It is inevitable that the particles will 
suffer adiabatic losses as the jets expand. We argue that 
synchrotron and inverse-Compton losses are negligible by 
comparison for electrons radiating at 8.4 GHz in 3C 31, as 
the jets have accurately p ower-law spectra with indices close 
to 0.5 ijLaing et alJ l2004'l and show no evi dence of spectral 
curva ture until much higher frequencies f Hardcastl e et alJ 
2002). It is therefore worthwhile to compare the observations 
with models in which the radiating particles are accelerated 
before entering the region of interest and then lose energy 
only by the adiabatic mechanism while the magnetic field 
is frozen into and convected passively with the flow, which 
is assumed to be laminar. Following conventional usage, we 
refer to such models as adiabatic. 

There are analytical solutions for the adiabatic evolu- 
tion of emissivity and surface brightness with jet radius and 
velocity if the magnetic field is exactly parallel or perpen- 
dicular to the flow d i rectio n. Models of this kind were first 
considered by iBurchl Jl978ft . who showed that the decrease 
of brightness with distance from the nucleus would be much 
steeper than he observed i n 3C 31 if the jets have constant 
velocity. iFanti et alJ lll982f) pointed out that the brightness 
would decline more slowly with distance in a dec elerating 
jet. This ide a underlies the turbulent jet models of [Bicknelll 
(1984, 1986). A rclativistic generalisation was developed by 
iBaum e t al.| 1)19 971) to mode l th e jets in 3C 264 and was ap- 
plied bv lFeretti et alJ (I1999T) and lBondi et alJ (|2000F) to other 
FRI jets. These treatments all assumed that there is no ve- 
locity gradient across the jet, and that the magnetic field 
is exactly parallel or perpendicular to the flow. While this 
approach is self-consistent, it cannot be applied if there is 
velocity shear in a direction perpendicular to any component 
of the field. 

In LB, we developed a numerical approach to modelling 
a relativistic jet with a more general magnetic field structure. 
This approach allowed us to use the variations of total in- 
tensity and linear polarization as independent constraints on 
the jet velocity field. We concluded that the magnetic struc- 
ture and velocity field in 3C 31 ar e indeed sign if icantl y more 
complex than those assumed by IBaum et alJ lll997f) . Nev- 
ertheless, we compared our models and the data for 3C 31 
with their analytical solutions. We showed that: 

(i) the adiabatic approximation is qualitatively inconsis- 
tent with the variations of brightness and polarization along 
the first «3 kpc of the jets, but 

(ii) further from the nucleus, the observed variations are 



closer to those expected if the adiabatic approximation 
holds. 

This comparison motivated us to develop a more general 
approach to modelling of adiabatic, relativistic jets, which 
we present in this paper. Our approach allows us to calculate 
how the brightness and polarization structure evolve along 
an adiabatic jet from prescribed initial conditions (specified 
as profiles across the jet), given the jet geometry and more 
complex magnetic structures and velocity fields of the type 
inferred for 3C 31 by LB. Detailed comparison of the new 
adiabatic models with the free models of LB, which fit the 
observations better, can then diagnose whether and where 
other physical processes, such as particle acceleration, may 
be significant in the jets. 

Sect ion [5] reviews our assumptions and the previous an- 
alytical solutions, and then outlines our calculation of syn- 
chrotron emission from adiabatic flows. Section [3] describes 
our approach to modelling of adiabatic jets; it briefly re- 
capitulates material from LB before discussing new aspects 
specific to the present study. Section[I]applies our adiabatic 
models to the outer regions of the jets in 3C 31 and shows 
that they can give a fair description of the VLA observa- 
tions of these regions. Section [S] confirms that the adiabatic 
models fail to describe the inner jet regions and critiques 
the adiabatic hypothesis in the light of this result; it also 
discusses the extent to which distributed particle injection 
can bring the adiabatic models into better agreement with 
the data. Section summarizes our conclusions. 



2 THE ADIABATIC APPROXIMATION 
2.1 Assumptions 

The jets are t aken to be adiabatic i n the sens e defined by 
IBurchl lll979l) . lMatthews fc Scheued<199d) and lBaum etall 
l|l997h . as follows: 

(i) The energies of the radiating particles change like 
those in an adiabatically expanding relativistic gas, i.e. 
oc I/ -1 / 3 , where V is the volume of a fluid element in its 
rest frame. 

(ii) There is no diffusion of particles. 

(iii) The particle momentum distribution remains 
isotropic (e.g. by resonant scattering off Alfven waves). 

(iv) The magnetic field behaves as if it is convected pas- 
sively with the flow. We take the velocity field to be a smooth 
function of position. Although the addition of a turbulent 
velocity component would make little difference to the cal- 
culation of the effects of relativistic aberration on the ap- 
pearance of the jet, there would be a major difference in the 
strength and structure of the magnetic field, which would 
be affected by shear and expansion, if not by dynamo action 
and reconnection. 

(v) Synchrotron and inverse-Compton losses are negligi- 
ble for electrons radiating at the wavelength of observation, 
so the energy and synchrotron frequency spectra are always 
power laws with indices —(2a + 1) and —a, respectively. 
Specifically, we take the number density of radiating elec- 
trons with energies between E and E + dE in the jet rest 
frame to be 

N(E)dE = nE- {2a+1) dE (1) 
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In addition, we assume, as in LB, that the regions of the 
jets to be modelled are intrinsically identical, antiparallel, 
axisymmetric, stationary flows. 

The magnetic field B is taken to have longitudinal, 
toroidal and radial components Bi, Bt and B r (all mea- 
sured in the jet rest frame). Synchrotron radiation is gener- 
ally anisotropic even in the rest frame of the emitting mate- 
rial. In this paper we therefore write the emissivity e(r)g(r), 
where e(r) would be the emissivity in total intensity for a 
magnetic field (Bf + B 2 + B 2 ) 1 ^ 2 perpendicular to the line 
of sight ({) denotes a spatial average). e(r) is the same for all 
three Stokes parameters, but g(r) depends on field geometry, 
and differs for /, Q and U [for total intensity, < gi(r) < 1 
and for linear polarization — po < gQ,u(r) < Po, where 
po = (3a + 3)/ (3a + 5) is the maximum degree of polar- 
ization for spectral index a]. 

2.2 Analytical approximations 

We now briefly recapitulate t he an alytical adiabatic rela- 
tions derived by IfBaum et alJ il997f) for axisymmetric, de- 
celerating, relativistic jets without velocity shear. Suppose 
that a jet has radius R, and that we can make the quasi-one- 
dimensional approximation: (a) that it is uniform in cross- 
section at a given distance from the nucleus and (b) that 
the velocity is unidirectional. The field components and the 
normalizing constant in the energy spectrum, n, vary as: 



Bi oc R~ z 

B t oc (R/3T)- 1 

B r oc (R/3T)- 1 

n oc (r/3ii 2 )- (1+2Q/3) 



(2) 
(3) 
(4) 
(5) 



where the velocity of the jet is j3c and T = (1 — /3 2 ) -1 ^ 2 . 

For a purely longitudinal field (B t — B r = 0), this leads 
to a variation of the rest-frame emissivity: 

-(2a+3)/3^-(10a+12)/3 



(IT*)" 



and for a perpendicular field (Bi = 0): 

6 OC (r/3) -(5-+6)/3 jR -(^+9)/3 



(6) 



(7) 



For conical jets, the radius is proportional to distance from 
the nucleus, r, so identical relations are obtained with r 
replacing R. 

In the absence of velocity shear, the adiabatic relations 
for the magnetic field can be combined to describe arbitrary 
initial conditions. The total field B is then: 

1/2 



B = 



{Bt + B'r 



(8) 



where Bi, Bt and B r are the initial field components, and 
F are the velocity and Lorentz factor, all defined where the 
jet radius is R = R. The emissivity is then: 



e oc {rpR 2 



-(l+So/3) nl+o 



(9) 



The next subsection develops a numerical approach capable 
of describing velocity shear. 



2.3 Adiabatic flows with arbitrary initial 
conditions and velocity fields 

The formalism needed to predict the total and linearly 
polarized synchrotron emission from an element of fluid 
in a non-relativistic a diabat ic flow was first developed by 
iMatthews fc Scheuerl jl99fj) . Their approach was to fol- 
low an element of fluid containing relativistic particles and 
an initially isotropic, disordered magnetic field through a 
model flow. They included the effects of synchrotron and 
invers e -Com pton energy losses, as well as adiabatic effects. 
lLaind (l2002t> developed a simplified approach for the case 
where synchrotron and inverse-Compton losses are negligi- 
ble, which we follow here. The reasons for taki ng the field t o 
be disordered on small scales are dis cussed bv lLainJ il98lf) . 
Bcgclma n. Blandford fc Reea 111984) and LB. 

An element of fluid is taken to be a unit cube (much 
smaller than the spatial scale of variations in the velocity 
field), with sides defined by orthonormal unit vectors a, b 
and c along natural axes of the flow in the fluid rest frame. 
The cube contains an isotropic disordered field and relativis- 
tic particles with the energy spectrum of equation Q and 
normalizing constant n. The cube moves with the flow, de- 
forming into a parallelepiped with sides defined by the vec- 
tors a, b and c. Application of flux conservation shows that 
an element of field which is initially B = B x a + B y h + B z c 
becomes: 



B 



B x a + B y h + B z c 
V 



(10) 



where V = a • b x c is the volume of the parallelepiped eval- 
uated in the fluid rest frame. The energy spectrum constant 
n evolves according to: 



n = nV 



-(l + 2a/3) 



(11) 



We assume that the flow is axisymmetric, with decel- 
eration and a transverse velocity gradient. We define a sec- 
ond coordinate system (x, y, z), again in the observed frame, 
with z along the jet axis, x perpendicular to it in a plane 
containing the line of sight and y in the plane of the sky. 
We parameterize the streamline by an index s, which varies 
from at the inner edge of a flow component to 1 at its outer 
edge, as in LB. Without loss of generality, we can consider 
flow in the x, z plane, where the distance of a streamline 
from the jet axis is x(z, s). 

The geometry and evolution of the vectors a, b and c are 
sketched in Fig. As an element of fluid moves outwards, c 
remains parallel to the streamline and its magnitude is deter- 
mined entirely by the flow velocity, b is unaffected by shear 
and is orthogonal to the flow direction, so its magnitude is 
proportional to the distance of the streamline from the jet 
axis, a (initially radial) is the only one of the three vectors 
affected by velocity shear: it is a function of the jet radius 
and the accumulated path length difference along adjacent 
streamlines (Fig. 0. Four quantities are therefore needed 
to describe the shape of the parallelepiped: three expansion 
factors and a shear term. For a streamline in the xz plane, 
the vectors are: 



dx/ds / i + x ,2 \ 1/2 r 
dx/ds I i + x' 2 J f /n 



b = 



x 

— m 

x 



(12) 
(13) 
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r/3 



and the volume is 

FB x dx/ds 
V = a ■ b x c = += 



Ffi x dx/ds \ l+x 



l+x'- 



1/2 



(14) 



(15) 



1, m and n are (ortho-normal) unit vectors along the local 
radial, toroidal and longitudinal directions, respectively and 
primes denote differentiation with respect to z for a given 
streamline. Barred quantities are evaluated at the starting 
surface. The factor of F/F in equations 1121 and 1141 accounts 
for Lorentz contraction along the flow direction. 

/ is the shear term, which does not affect the volume. It 
will usually be < if the velocity decreases outwards from 
the jet axis. We evaluate it by considering two elements of 
fluid which leave the reference surface at time t = on 
adjacent streamlines with indices s and s + As. After a time 
t, we have 



ct 



1 



' As) [l + x'(z,s) 2 ]^ 2 dz 



Z(s+As) [l + x'(z,s + As) 2 ]^ 2 dz 



z(s + As) 



(3{z, S + As) 



(16) 

where dl is an element of path along the streamline and z(s) 
and z(s) are the z-coordinates of the fluid element at the 
reference surface and after time t, respectively, We express 
the integral for streamline s + As as the sum of the integral 
for streamline s and a set of terms expanded to first order 
in As and Az — z(s + As) — z(s). The reference surface is 
spherical, and is always set in a part of the jet where the 
streamlines are straight, so a streamline crosses the surface 
at z = z with 



x(z, s) = r cos[C + (| — C) s l 



(17) 



(see Section I3.2H . This allows us to evaluate Az/As using 
the relation 



Az 
As" 



l+x'- 



($-C)ftan[C + (g-Qs] 

$ 

f3x'dx'/ds - (1 + x' 2 )df3/ds 
f3 2 (l + x' 2 ) 1 / 2 



dz 



(18) 

The difference in path length along the streamline, Al, is 
given by 



— = (1+3! 

As 



/2\l/2 



Az 
As" 



dx 



ds 1 + x' : 



(19) 



Finally, the shear term / is the component of a along the 
streamline normalized by the magnitude of the vectors at 
the reference surface: 

'"^ 

/ is negative if the velocity decreases with increasing s. 

In general, calculation of the shear term requires a nu- 
merical integration, but for regions of the jet where the 



streamlines are straight and the variation of the velocity 
along a streamline is a simple function, it can be done ana- 
lytically. Note also that the shear term is non-zero even for 
a velocity independent of s if the streamlines are curved. 

For flow radially outwards from the nucleus, the three 
vectors and the volume element can be written very simply 
in terms of the distance from the nucleus, r — (x 2 + y 2 + 

J2W2. 



V = 



r. r 
F a+ f /c 

r 

V M 

Ff3 

FJ3f_ 

FPr 2 



(21) 
(22) 
(23) 

(24) 



In the absence of sh ear (/ = 0) these a re equivalent to the 
relations derived bv lBaum et 

aD <ll997fl and given m equa- 
tion s HE] 

iLaind (120021) gives the variation of the field components 
perpendicular to the line of sight in terms of the direction 
cosines of the vectors a, b and c with respect to a fixed 
coordinate system (X, Y, Z) with Z along the line of sight. 
Together with n, these determine the synchrotron emissivity 
in Stokes 7, Q and U. 

There are two additional complications: 

(i) The direction cosines must be evaluated in the rest 
frame of the fluid (exactly as for the free models of LB) . 

(ii) We need to set initial conditions on a reference surface 
in the jet, rather than starting with an isotropic magnetic 
field. We do this by adjusting the relative magnitudes of the 
vectors a, b and c in such a way as to produce the desired 
ratios between the three field components whilst leaving the 
volume (and therefore the rms total field and the emissivity 
function, e) unchanged. 



3 OBSERVATIONS AND MODELLING 
METHODS 

3.1 Observations 

The deep, high-resolution VLA observations with which we 
compare our models are described in LB. We take the Hub- 
ble constant to be Hp = 70 km s" 1 Mpc" 1 . At the redshift of 
3C31/0.0169:ISmith et alJ[200d|Huchra. Voeelev fc Gelled 

~w7~ 



19991 iDe Vaucouleurs et al.lll99 



the linear scale is then 
0.34 kpc/arcsec. We fit our models to images at resolutions 
of 0.75 and 0.25arcsec, covering the inner ±28arcsec of the 
jets. Fig. |5| shows the emission from the jets of 3C31 with 
the area we model outlined. 



3.2 Geometry 

In LB, we showed that the jets in 3C31 could be divided 
into three regions according to the shapes of their outer 
isophotes. As observed (i.e. projected on the plane of the 
sky), these are: 

(i) Inner (0 - 2.5 arcsec): a cone, centred on the nucleus, 
with a projected half-opening angle of 8.5 degrees. 
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/ . a \ 






/s + As \ 
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a \ 


" s 




C 








Figure 1. Sketches of the evolution of the vectors a, b and c de- 
fined in the text (not to scale). Panel (a) shows two neighbouring 
streamlines with indices s and s + As (defined in Section 12.31 . in 
a plane containing the jet axis in a part of the jet where the flow 
direction changes with distance from the nucleus, b and b are 
perpendicular to the plane of the diagram, and are not shown. 
Note that c remains parallel to the flow despite the change in 
direction. Panel (b) shows a projection of the flow in three di- 
mensions viewed with the jet axis at 45° to the line of sight in a 
region where the streamlines are straight. The full line is a repre- 
sentative streamline. The ellipses represent the positions of fluid 
elements with streamline indices s and s + As. The lower pair are 
at the surface where the initial conditions are set and the upper 
pair represent the positions reached by the fluid elements after a 
given time. The vector a', shown dashed, is the value of a in the 
absence of shear. 




Figure 2. A grey-scale of the total intensity from the jets of 
3C31 at 8.4 GHz (LB). The resolution is 0.25arcsec FWHM and 
the area we model is indicated by the box. 



(ii) Flaring (2.5 - 8.3 arcsec): a region in which the jet 
initially expands much more rapidly and then recollimates. 

(iii) Outer (8.3 - 28.3 arcsec): a second region of conical 
expansion, also centred on the nucleus, but with a projected 
half-opening angle of 16.75 degrees. 

We showed that these regions also have distinct kinematic 
properties. 

We use the same descriptions of jet geometry and veloc- 
ity as in LB, where we investigated two different transverse 
velocity structures. In the first {spine/shear layer - SSL) 
a central fast spine with no transverse variation of velocity 
is surrounded by a slower shear layer with a linear velocity 
gradient. In the second (Gaussian), there is no distinct spine 
component, and the jet consists entirely of a shear layer with 
a truncated Gaussian velocity law. The streamline index, s 
(defined in table 3 of LB) varies from for the streamline 
closest to the axis in the spine or shear layer to 1 for the 
furthest streamline. The angles f and | used in Section [2.31 
are the opening angles of the spine and shear layer, respec- 
tively, at the reference surface for an adiabatic model. 
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Table 1. The functional variation of velocity along the model jets. 



Quantity Free parameters 

ri ro Tf Other 



Functional dependences 
Flaring Outer 



Mp) 

v{p) 



On-axis velocities 
/3i Po ft H bo + hpH- 1 



+ b 2 p H c exp(-C!p) 



<-'(> 



Fractional edge velocities 

(p-r 1 )(v -v 1 ) 

r —ri 



'■•It 




(a) 



In order to describe velocity variations along a stream- 
line, we use a coordinate p, defined as: 



(b) 


^^^^ ^^^^ " " 










r. --^ 
Inner 


~ ~ - 

Flaring / 
5 arcsec 





Figure 3. Geometry of the spine/shear-layer model, showing the 
inner, flaring and outer regions in the plane containing the jet 
axis. The thick full curves represent the edge of the jet, the bound- 
aries between regions are represented by thin full curves and the 
s = 0.5 streamlines for the spine and shear layer are drawn as 
dashed curves, (a) The entire modelled region; (b) the base of the 
jet on a larger scale, showing the boundary surfaces at distances 
of n and ro from the nucleus. The Gaussian model is essentially 
the same, but with the spine component removed. 



The angle between the jet axis and the line of sight is 
taken to be 9 and we use the xyz coordinate system defined 
earlier. In the inner and outer regions the streamlines are 
straight. In the flaring region we interpolate using a cubic 
polynomial in such a way that x(z,s) and its first spatial 
derivative x'(z, s) are continuous at the boundaries between 
regions, which are spherical, centred on the nucleus with 
radii n and ro. The assumed geometry is sketched in Fig. 01 



p — r (inner region 
p = n + (ro — ri) 



ri cos 4>i(s) 



ro cos <t> (s) — ri cos (pi(s) 
(flaring region) 
r (outer region) 



where the streamline makes angles <f>i and (j> with the axis 
in the inner and outer regions, respectively, p is monotonic 
along any streamline and varies smoothly from ri to ro 
through the flaring region. This allows us to match on to 
simple velocity profiles which depend only on r in the inner 
and outer regions. 

The functions defining the edge of the jet are con- 
strained to match the observed outer isophotes and are fixed 
in a coordinate system projected on the sky. Their values in 
the jet coordinate system then depend only on the angle to 
the line of sight. 



3.3 Velocity field 

The velocity field is taken to be a separable function 
(3(p,s) — f3 p (p)/3 s (s) with p s (0) = 1. The variation along 
a streamline, /3 p (p) is given in Table Q and is exactly as 
in LB, omitting the inner region, which we do not discuss 
quantitatively in this paper. It is defined by the index H, 
together with velocities at three locations in the jet: n, ro 
and an arbitrary fiducial distance r/. These distances are 
fixed by their projections on the plane of the sky, which are 
2.5 arcsec, 8.2 arcsec and 22.4 arcsec, respectively. The form 
of /3 p (p) was chosen by LB (Section 3.4) to permit fitting 
the sidedness-ratio profile observed in 3C31, which requires 
the velocity to remain constant through most of the flaring 
region but then to drop abruptly close to the outer bound- 
ary, at a rate determined by the index H. The velocity then 
falls smoothly and slowly through the outer region. 

/3 a (s) = 1 in the spine for SSL models, dropping lin- 
early with s from 1 at the spine/shear layer interface to a 
minimum value at the edge of the jet. For Gaussian models, 
@s(s) is a truncated Gaussian function. In all models, the 
minimum fractional velocity, v(p), is allowed to vary along 
the jet (Table [TJ. 
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Table 2. Summary of the functional variations of velocity, emis- 
sivity, and field component ratios across the jet. 



Location Functional variation 



Transverse velocity profile (varies along jet) 
Spine /3 s (s) = l 

Shear layer SSL /3 s (s) = 1 + [v(p) - l]s 

Shear layer Gaussian /3s(s) = exp[— s 2 mi)(p)] 

Emissivity profile at p = r = f 
Spine e(f , s) = E 

Shear layer SSL e(f , s) = 1 + (e - l)s 

Shear layer Gaussian e(f, s) = exp(s 2 lne) 

Radial/toroidal field ratio at p = r = f 
Spine j(f, s) = J spinc 

Shear layer j(f , s) = j ccn + s p (j out - j ccn ) 

Longitudinal/toroidal field ratio at p = r = f 
Spine k(f, s) = fc sp i ne 

Shear layer fc(f, s) = fc ccn + s q (k out - Scon) 



(i) Determine coordinates in a frame fixed in the jet, in 
particular the distance coordinate p and the streamline in- 
dex s, numerically if necessary. 

(ii) Evaluate the velocity at that point, together with the 
angle between the flow direction and the line of sight ip. 
Derive the Doppler factor D = [F(l — /3cost/j)] -1 and hence 
the rotation due to aberration (sin?// = Dsini/>, where tp' is 
measured in the rest frame of the jet material). 

(iii) Look up the initial values for the emissivity and field- 
ordering parameters on the streamline. 

(iv) Calculate the evolution of the particle density and 
field components along the streamline using equations II II 
nillTHl and tTglOrjl 

(v) Evaluate the emissivity function e(p, s) oc nB +a . 

(vi) Evaluate the position angle of polarization, and the 
rms field components along the major and minor axes of 
the probability de nsity functio n of the field projected on the 
plane of the sky jLaind lFoC^. Multiply by e(p, s)D 2+a , to 
scale the emissivity and account for Doppler beaming. 

(vii) Derive the to tal and polar ized emissivities using the 
expressions given bv iLainsI (l2002f> and convert to observed 
Stokes Q and U. 

Fitting to observed images and optimization of the mod- 
els are done by x 2 minimization as described in LB. 



3.4 Initial conditions for emissivity and field 
ordering 

For a given velocity field in the jet, the adiabatic models 
require the radiating particle density and magnetic field 
components to evolve self-consistently with the velocity. 
These models can therefore be specified completely by set- 
ting their initial values at one point on any given stream- 
line, most straightforwardly at a surface of constant dis- 
tance from the nucleus p = r = f. We need to define the 
emissivity variation e(f , s) and two field-component ratios: 
j(f,s) = (B 2 ) 1/2 /(B 2 ) 1/2 (radial/ toroidal) and k(f,s) = 
(Bf) 1/2 /{B 2 ) 1/2 (longitudinal/toroidal). The number of free 
parameters needed to specify an adiabatic model is much 
smaller than that for the equivalent free model, in which the 
emissivity and field ordering are allowed to vary separately 
as smooth functions of position. 

We have made two sets of adiabatic models, attempting 
to fit different regions of the jet, as follows: 

(i) the outer region alone, with initial conditions set at its 
boundary with the flaring region (f = ro); 

(ii) the flaring and outer regions, with initial conditions 
set at the boundary between the inner and flaring regions 
(f = ri). 

The functional forms for the initial transverse variations of 
emissivity, e(f, s) and the field component ratios j(f, s) and 
fc(f, s) are given in Tabled 

3.5 Model integration, fitting and optimization 

The integration through the model jets to determine the I, Q 
and U brightness distributions is identical to that described 
by LB, except for the determination of emissivity and field 
ordering at a given point in the jet, the main steps in which 
are: 



4 FITS TO THE OUTER REGION ALONE 

We know ( LB, Fig. 20) that th e analytical adiabatic approx- 
imation of Ba um et alJ ((1997) is within a factor of two of 
the variation of emissivity along the outer jet if the velocity 
profile is as estimated in our best-fitting free models. In this 
section, we therefore fit full adiabatic models to the outer re- 
gion alone, setting the initial conditions at its boundary with 
the flaring region and and computing the \ 2 value for the fit 
at distances >10arcsec from the nucleus, where all lines of 
sight pass only through the outer region. We initially allow 
the angle to the line of sight, the velocity field and the initial 
emissivity and field- ordering parameters to vary. The results 
for the Gaussian and SSL adiabatic models are very similar, 
but the former always fit slightly better. As the Gaussian 
models also have fewer free parameters, we concentrate on 
them when comparing two classes of adiabatic models with 
the data, as follows: 

(i) A model having precisely the same angle to the line 
of sight, velocity field, initial conditions and flux scaling as 
the best-fitting free Gaussian model from LB, which it is 
therefore forced to match at the end of the flaring region. We 
henceforth refer to this as the parameter-matched adiabatic 
model. It has no free parameters. 

(ii) A fully-optimized adiabatic model. The model flux 
density is constrained to be the measured value for the outer 
region, but all other parameters are varied to achieve the 
best fit. 

For the best-fitting free model, there is a well-constrained so- 
lution for which the reduced x 2 m the outer region, Xicd ~ 
1.6 (LB). For the adiabatic models, the optimization routine 
failed to find a well-defined x 2 minimum: there is a broad 
range of solutions with very similar brightness and polariza- 
tion distributions, but different angles to the line of sight. 
We therefore fixed 9 at values separated by 5° in the range 
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Figure 4. Contours of total intensity and magnetic-field vectors at a resolution of 0.75 arcsec for the outer regions of the jet and 
counter-jet. The plots cover the range 10 - 27 arcsec on either side of the nucleus and are in pairs with the counter-jet on the left and the 
main jet on the right. The angular scale is indicated by the labelled bars at the bottom of the diagram. From the top: observations; free 
SSL model; parameter-matched adiabatic model; fully-optimized adiabatic model. The initial conditions for the adiabatic models are set 
at the boundary between the flaring and outer regions. The left-hand pairs of panels show total intensity contours at levels of —1, 1, 2, 4, 
8, 16, 32 X 20 /jjy/beam area. The right-hand pairs show vectors whose magnitudes are proportional to p and whose directions are those 
of the apparent magnetic field, superimposed on selected / contours. The scale of p is indicated by the bar in the bottom right-hand 
panel. 

20° < 9 < 75° and optimized the remaining parameters. not performed an error analysis because of the wide range 
The minimum Xrcd = 1-9 for the Gaussian model occurs of acceptable parameters, 
for ~ 55°, but there are solutions with Xrcd < 2.5 for 
25° < 9 < 70°. The best-fitting values are listed in Table 01 
along with those of the parameter-matched model. We have 
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Figure 5. Grey-scale representations of the degree of polariza- 
tion, p, in the outer jets. The grey-scale runs from p = to p = 0.7, 
as indicated by the labelled wedges. The right- and left-hand sets 
of panels show the main and counter-jets, respectively, and the 
angular scale is indicated by the labelled bar at the bottom of the 
diagram. From the top: observations; free SSL model; parameter- 
matched adiabatic model; fully-optimized adiabatic model. 
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Figure 6. Profiles of total intensity along the jet axes at a reso- 
lution of 0.75 arcsec FWHM. Full line: observations; short dashed 
line: free SSL model; dotted line: parameter-matched adiabatic 
model; long-dashed line: optimized adiabatic model. The initial 
conditions for the adiabatic models are set at the end of the flaring 
region, so their profiles are plotted only beyond 10 arcsec from the 
nucleus, where the line of sight is entirely within the outer region. 
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Table 3. Parameters for adiabatic models of the outer region. 
The initial values are denned at the boundary between the outer 
and flaring regions, f = tq. 



Parameter 



Parameter- Fully- 
matched optimized 



Geometry 
Angle to line of sight 9 

Velocity field 

On-axis velocity /3o 
On-axis velocity 0f 
Fractional edge velocity vg 



51.4 



0.54 
0.27 
0.63 



55.0 



0.76 
0.38 
0.26 



(iii) Both adiabatic models predict too high a degree of 
polarization in the transverse-field region on the axis of the 
main jet, especially within ~10arcsec of the start of the 
outer region (Fig. 01. 

(iv) Conversely, both adiabatic models give slightly too 
low a degree of polarization on the axis of the counter-jet 
(Figs. El and EJ. 

(v) The degree of polarization at the edges of both jets 
is overestimated, again within «10arcsec of the start of the 
outer region. (Figs. 2] and [KJ 

(vi) The field-vector directions for the fully-optimized 
adiabatic model are close to those observed, but those for 
the parameter-matched model are incorrect at the edge of 
the jet (where they should be parallel to the surface; Fig.0J. 



Emissivity profile 

Fractional edge emissivity e 0.26 0.26 

Radial/toroidal field ratios 

On-axis j ce n 0.0 0.53 

Edge j out 0.92 0.0 

Index p 0.41 1.09 

Longitudinal/toroidal field ratios 

On-axis k ccn 0.82 2.53 

Edge fcout 0.82 0.0 

Index q 0.0 0.55 

Goodness of fit 

X? od 8.64 1.89 



Fig. 2] compares the data and the free and adiabatic 
models by showing contours of total intensity for the outer 
regions of the jets, the best-fitting free SSL model and the 
two Gaussian adiabatic models. It also plots vectors whose 
magnitudes are proportional to the degree of polarization, 
p, and whose directions are those of the apparent magnetic 
field. Fig. shows p for the same regions in a grey-scale 
representation and Fig |5] shows longitudinal profiles of total 
intensity. 

The main features evident from this comparison are as 
follows: 

(i) The total intensity predicted by the parameter- 
matched adiabatic model falls off rapidly close to the bound- 
ary between the flaring and outer regions (where it is forced 
to match the free model) and thereafter is much lower than 
the observed values. This is most clearly shown by the 
longitudinal profile in Fig. |S] and is reflected in the high 
X?cd=8.6. 

(ii) The total-intensity distribution for the fully- 
optimized adiabatic model is in much better agreement 
with the observations. Some of the improvement in \ 
results from the requirement to fit the total flux density in 
the outer region rather than to match the free model exactly 
at the start of the region. The model parameters can then 
be optimized to give a good fit between 11 and 25arcsec 
from the nucleus, at the expense of an error between 10 
and llarcsec, where the initial decrease of brightness with 
distance from the nucleus is slightly too rapid, and at 
distances > 25 arcsec, where the model intensity is too 
high (Figs [I] and [§J. The total intensity of the counter-jet is 
well fitted by this model. 



The problems encountered in fitting the polarization of 
the outer region with these laminar-flow adiabatic models 
are fundamental and do not depend on the choice of 9 or 
the velocity field. The free models fitted by LB introduce 
a significant radial magnetic field component, B r , at the 
edge of the jet at the start of the outer region to produce 
the low degree of polarization p there. This field component 
also reduces the values of p on-axis. The very high degrees of 
polarization observed at the jet edges further from the nu- 
cleus require B r again to become small compared with the 
toroidal and longitudinal components Bt and Bi there. This 
cannot be achieved in our adiabatic models, for which the 
ratio B r /Bt is fixed by our choice of an axisymmetric veloc- 
ity field with straight streamlines for the outer region (the 
coefficients of a and b in equations 1211 and 1221 are identical). 
More generally, B r /B t alters very slowly unless the rate of 
change of expansion is large: this is a direct consequence of 
flux-freezing in a laminar velocity field (equations 1121 and 
1131 . In order to change the ratio, it is necessary to introduce 
a component of velocity across the existing model stream- 
lines. This component could result from turbulence, for ex- 
ample as a result of entrainment of the external medium. 
The reason why B r disappears in the outer parts of the jets 
remains unclear. 

As a result of the small range of jet/counter-jet sided- 
ness ratio in the outer region, the fully-optimized adiabatic 
models are poorly constrained. This results in a degeneracy 
between the velocity and angle to the line of sight: the opti- 
mized values of the fiducial velocities fio and /3f range from 
0.33 and 0.16 for 9 = 25° to 0.87 and 0.43 for 9 = 70°. For the 
best-fitting model, the velocities are 0.71 and 0.34 (Tabled . 
These are significantly larger than the velocities deduced for 
the free model (0.54 and 0.27) and imply faster deceleration. 
In addition, the fractional velocity at the edge of the jet is 
lower (0.26, compared with 0.63). The combination of these 
two differences leads to much higher shear at the edges of 
the jet for the optimized model. The changes result from a 
partially-successful attempt to fit the slow brightness decline 
at the start of the outer region and the apparent magnetic- 
field direction at the edge of the jet. In effect, the coupling 
that the adiabatic models require between the variations of 
the field geometry and the emissivity forces them to steeper 
velocity gradients in the outer region than those in the best- 
fitting free model. We have also investigated the effects of 
changing the functional form of the velocity law as well as 
the parameters. It is possible to improve the fit to the total 
intensity by making the velocity decrease more rapidly at 
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the start of the outer region and more slowly at large dis- 
tances, but this always gives a worse fit to the polarization, 
and we have been unable to improve the overall \ 2 - 

The parameter-matched model effectively incorporates 
information gained from free-model fits to the inner jet re- 
gions and should therefore provide a much more realistic 
description of the velocity field, but its predicted bright- 
ness distribution initially declines too rapidly with distance 
from the nucleus (Figs 0] and EJ. Particle injection and/or 
field amplification might occur close to the boundary with 
the flaring region, thereby slowing the emissivity decline (cf. 
Section 



5 FITS THAT INCLUDE THE FLARING 
REGION 

5.1 The inner region 

As noted by LB (Section 5.4), the conical inner region (the 
first 2.5 arcsec of the jet) shows no evidence for decelera- 
tion, but neither does it exhibit the extremely rapid bright- 
ness fall-off characteristic of adiabatic expansion at constant 
speed. The low sidedness ratio observed in this region also 
led us to suggest that its emission comes mostly from a slow 
surface layer which does not persist to larger distances. We 
have insufficient resolution to build or verify a model of the 
transverse structure of this region, so we do not consider it 
further in this paper. 

5.2 Models with initial conditions set at the 
boundary between the inner and flaring 
regions 

LB showed t hat th e analytical adiabatic approximation of 
iBaum et al.l (Il997f) fails completely for the flaring region. 
We now examine whether this conclusion would be modified 
by including more realistic initial conditions and the effects 
of velocity shear by modelling the flaring and outer regions 
of 3C 3f , setting the initial emissivity and field- ordering pro- 
files at the boundary between the inner and flaring regions 
(f = n). As in Section 0] we find that the results for Gaus- 
sian and SSL velocity profiles are very similar, but that the 
former fit better, as well as having fewer free parameters. We 
therefore again show only the Gaussian-velocity case, and 
compare fully-optimized models with a parameter-matched 
model whose initial conditions, velocity field, angle to the 
line of sight and flux scaling are identical to those of the 
best-fitting free model. 

The best fits for the optimized models are much poorer 
(Xred ~ 5) than those for the free models (x? G d — 1-7 _ L8), 
and also require extremely high velocities at the start of 
the flaring region. As for the outer region alone (Section 
solutions of comparable quality can be found over a wide 
range of angles to the line of sight (30° < 8 < 55°). We 
show the best-fitting example, which again has 6 = 55°. 

Images of total intensity and linear polarization for 
the parameter-matched and fully-optimized models are com- 
pared with the observations of 3C 31 and the best-fitting free 
SSL model in FigsfTl- llUl and the parameters for both mod- 
els are listed in Tabled The main points of interest are as 
follows. 



(i) The parameter-matched model fails completely to fit 
the brightness distribution: its initial brightness fall-off is 
far too steep on both sides of the nucleus. This confirms our 
conclusion from the analytical approximation. 

(ii) The polarization distribution predicted by the 
parameter- matched adiabatic model is also incorrect, the ap- 
parent magnetic field being almost perpendicular to the jets 
at their edges, rather than parallel as observed. 

(iii) The fully-optimized model fits the total intensity 
from the main jet and the outer counter-jet fairly well 
(Figs 171 l9land rT"(H . but seriously underestimates the emission 
from the counter-jet close to the beginning of the flaring re- 
gion, where the predicted jet/counter-jet sidedness ratio is 
~50, compared with the observed maximum of 13 for the 
entire region. 

(iv) The fully-optimized model shows a qualitatively cor- 
rect polarization distribution, with transverse apparent field 
on the axis of both jets and longitudinal field at the edges. 
It predicts a higher degree of polarization on-axis in the 
main jet than in the counter-jet (opposite to the observed 
difference) and also fails to reproduce the regions of low po- 
larization at the edges of both jets and across the whole of 
the main jet in the flaring region (Figs Q and |5J. 

The parameter-matched model fails because there is in- 
sufficient deceleration to counteract adiabatic losses result- 
ing from the rapid expansion of the jets in the flaring region, 
even when field amplification by velocity shear is taken into 
account. There is also insufficient shear to counteract the 
effects of expansion at the edges of the jets, leading to a 
transverse apparent field there. The fully-optimized model 
fits much better, but has two main problems. First, it re- 
quires an extremely rapid deceleration in the flaring region: 
the initial on-axis velocity f3i is in the range 0.93 - 0.99 for 
any 9 and the corresponding velocity at the edge of the jet 
is > 0.7, so high sidedness ratios are inevitable. Second, as 
in the outer region, the polarization data require far larger 
a change in the ratio B r /Bt than is allowed in our axisym- 
metric, laminar adiabatic model: a significant radial field 
component must first be generated at the edge of the flaring 
region and then destroyed further from the nucleus. 

The extremely high sidedness ratios required by the op- 
timized adiabatic models are inconsistent with observations. 
To help adiabatic models fit the data using the more real- 
istic velocity field of the free models, additional emissivity 
must be introduced in the flaring region, as we now discuss. 

5.3 Particle injection 

ft is clear from the earlier discussion that processes other 
than adiabatic evolution in an axisymmetric, laminar veloc- 
ity field must be important in the flaring region of 3C 31 and 
may also have some effect at the start of the outer region. A 
process that locally increases the emissivity is clearly needed 
to compensate adiabatic expansion losses in the flaring re- 
gion. Our modelling approach lets us estimate the rest-frame 
emissivity and the field-component ratios, but it cannot dis- 
entangle the roles of particles and magnetic field in the ab- 
sence of any constraint from inverse-Compton emission. The 
observation of X-ray emission from the inner and flaring re- 
gions shows that fresh relativistic particles must be injected 
there, as the synchrotron lifetimes of electrons radiating at 
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Figure 7. Contours of total intensity and magnetic-field vectors at a resolution of 0.75 arcsec. The plots cover 27arcsec on either side 
of the nucleus and the angular scale is indicated by the labelled bars. From the top: observations; free SSL model; parameter-matched 
adiabatic model; fully-optimized adiabatic model; adiabatic model with distributed particle injection. The initial conditions for both 
adiabatic models are set at the boundary between the inner and flaring regions and the inner region is not modelled. The left-hand panels 
show total intensity contours at levels of —1, 1, 2, 4, 8, 16, 32 X 20 pjy/beam area. The right-hand pairs show vectors whose magnitudes 
are proportional to p and whose directions are those of the apparent magnetic field, superimposed on selected / contours. The scale of 
p is indicated by the bar in the bottom right-hand panel. No emission is calculated in the inner region (±2.5 arcsec) for the adiabatic 
models. 



these frequencies are ~10's of years l|Hardcastle et afl 2002 'l. 
We therefore investigate a simple, but self-consistent physi- 
cal model in which new (e.g. reaccelerated) particles with an 
energy distribution An(E)dE oc £T (2a+1) <iE are added in 
a distributed fashion over the flaring region and then evolve 
adiabatically. We again fix the angle to the line of sight and 
velocity field at the best-fitting free model values, and as- 
sume that the magnetic field is frozen into the flow. We 
then optimize the initial conditions and particle injection 
function. We parameterize the particle injection as the ratio 
of the additional emissivity at position (p, s) in the flaring 
region to its value on the same streamline at the start of the 
region, using the empirical relation; 



Mft£) 
e(ri,s) 



= A 



p — r\ 

ro — n 



for n < p < pinj (the emissivity function e is defined in 
Section fe.lll . This means that the transverse injection profile 
is a constant Gaussian function. The parameters An, Q and 
pinj are optimized, with pi n j restricted to n < pi n j < ro. 

The fit to the total-intensity distribution is generally 
good (Figs 171 l9l and [TUll . but the jet/counter-jet sidedness 
ratio close to the boundary between the inner and flar- 
ing regions is still significantly higher than is observed, de- 
spite the more modest velocity there. The reason is that 
the field is almost purely lon gitudinal, so th e ratio is ~ 
[(l + /3cos0)/(l-/3cos6O] 3+2a feegelmarll995t rather than 
the usual « [(1 + /3cos#)/(l — [3 cos 6)] 2+a appropriate for 
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Figure 8. Grey-scale representations of the degree of polariza- 
tion. The grey-scale runs from p = to p = 0.7, as indicated by 
the labelled wedge. The panels extend ±27arcsec from the nu- 
cleus and the angular scale is indicated by the labelled bar. From 
the top: observations; free SSL model; parameter-matched adi- 
abatic model; fully-optimized adiabatic model; adiabatic model 
with distributed particle injection. The initial conditions for the 
adiabatic models are set at the boundary between the inner and 
flaring regions. No emission is calculated for the adiabatic models 
within 2.5arcsec of the nucleus. 



a field which is closer to isotropy. Given that the assumed 
magnetic-field evolution is unchanged, it is not surprising 
that the deficiencies in the polarization fits noted earlier are 
still present (Figs|7|and|HJ and that an unrealistic field struc- 
ture is required at the start of the flaring region. We do not, 
therefore, regard the high sidedness ratio as a significant 
problem for the model. The fit to the total intensity in the 
outer region is very similar to that of the parameter-matched 
model described in Section [I] as expected, so the need for 
some additional emissivity there remains. The overall fit is 
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Figure 9. Profiles of total intensity along the jet axes at a reso- 
lution of 0.75 arcsec FWHM. Full line: observations; short dashes: 
free SSL model; dots: parameter-matched adiabatic model; long 
dashes: fully-optimized adiabatic model; dash-dots: adiabatic 
model with distributed particle injection. The initial conditions 
for the adiabatic models are set at the boundary between the in- 
ner and flaring regions, so their profiles are plotted only for the 
outer and flaring regions (> 2. 5 arcsec from the nucleus). 



slightly worse than for the fully-optimized model (x? c d = 4.9 
compared with 4.7), but with 10 free parameters instead 
of 14. The optimized parameters for the particle-injection 
model are given in Table |3] Key implications are as follows: 

(i) Particle injection is required out to 6.4 arcsec (2.2 kpc) 
from the nucleus, i.e. only over the first half of the flaring 
region. In the second half of the region, the jet decelerates 
rapidly, thereby slowing the brightness decline predicted by 
adiabatic models. Further particle injection there would pro- 
duce too slow a decline in brightness. 

(ii) The injection rate decreases rapidly with distance 
from the nucleus within that region. 

(iii) The emissivity everywhere in the flaring region is 
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Figure 10. Contours of total intensity at a resolution of 
0.25 arcsec FWHM. The panels cover ±10 arcsec from the nucleus 
and the angular scale is indicated. The contour levels are —1, 1, 2, 
4, 8, 16, 32 X 30 ^tjy/beam area. From the top: observations; free 
SSL model; parameter-matched adiabatic model; fully-optimized 
adiabatic model; adiabatic model with distributed particle injec- 
tion. As the adiabatic models have initial conditions set at the 
boundary between the inner and flaring regions, no emission is 
calculated for the inner jet region (<2. 5 arcsec from the nucleus). 



Table 4. Fitted parameters for models with initial conditions at 
the boundary between the inner and flaring regions (r = n). Note 
that 9 and the velocity field are fixed for the particle-injection 
model; its remaining parameters are optimized. 



Quantity 



Parameter- Fully- Particle 
matched optimized injection 



Angle to line of sight 



9 (degrees) 


51.4 


55.0 


51.4 




Velocity field 






On-axis velocities 










0.76 


0.99 


0.76 


Pa 


0.54 


0.77 


0.54 


Pf 


0.27 


0.40 


0.27 


Edge velocities 








VI 


0.97 


0.93 


0.97 


Vo 


0.63 


0.30 


0.63 


velocity exponent H 


8.82 


2.00 


8.82 




Emissivity profile 






Edge emissivity e 


0.37 


0.05 


0.18 


Radial/toroidal field ratios 






On-axis j ce n 


0.0 


0.76 


0.81 


Edge j ou t 


0.78 


0.42 


0.42 


Index p 


0.41 


1.50 


0.62 


Longitudinal/toroidal field ratios 




On-axis fc ce n 


1.17 


4.97 


3.01 


Edge fc ut 


1.17 


2.65 


3.01 


Index q 


0.0 


3.96 


0.0 



Particle injection parameters 



Cut-off pi n j (kpc) 
Exponent Q 
Normalization Aq 



X 2 , 

A red 



Goodness of fit 
30.4 



4.65 



2.16 
-4.71 
84.5 



1.92 



dominated by particles injected locally rather than those 
originating from the inner region. 

Given the short synchrotron lifetimes of the electrons re- 
sponsible for the high-energy emission, X-rays are likely to 
be produced close to the injection sites and the X-ray emis- 
sivity might therefore be expected to be roughly propor- 
tional to the injection function, A, rather than to the radio 
emissivity. A compariso n between the observed X-ray and 
radio emission profiles (IHardcastle et al.l I2002L LB) and a 
model profile derived on this assumption (Fig. 111(1 supports 
this idea. The brightest X-ray emission comes from the in- 
ner 2.5 arcsec of the jet, which we have not modelled, but 
for which non-adiabatic evolution is clearly required (Sec- 
tion l5.1|l . In the first half of the flaring region, where we infer 
the need for significant particle injection, there is indeed ev- 
idence for more X-ray emission than in the outer half of the 
region, where injection is not required. The X-ray data in 
the flaring region are, however, too noisy to be sure whether 
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the ratio of X-ray to radio emission changes within this re- 
gion as our particle-injection model predicts. We also note 
that high-resolution studies of X-ray and radio emission from 
jets have shown that both may be highly inhomogeneous , 
with a complex mutu al relationship llHardcastle et af]l2003t 
IWilson fc Yanell2002ft . 

If particles are indeed injected without much modify- 
ing the magnetic field, the ratio of particle to field energy is 
likely to increase (by an amount which we cannot estimate 
without a proper description of the injection mechanism). 
We would not, therefore, expect equipartition between field 
and particle energy everywhere in the flaring region, even if 
it holds in some locations. The region where we infer sig- 
nificant particle injection is also where the expansion rate 
of the j et is increasing (Fig. and is where iLaine fc Bridle! 
(2002b) inferred a significant over-pressure (Fig. 111( 3 ). The 
internal pressure significantly exceeds not only the external 
gas pressure, but also the synchrotron minimum pressure, 
so quite large deviations from equipartition could occur be- 
tween 2.5 and «5 arcsec from the nucleus. In any case, we 
know that the magnetic field component ratios change in a 
way which is inconsistent with flux freezing in a laminar ve- 
locity field, so it is very likely that amplification of the field 
also occurs. 

All of these results are consistent with injection (or reac- 
celeration) of radiating particles over a wide energy range in 
the first half of the flaring region. The X-ray-emitting parti- 
cles have extremely short synchrotron lifetimes and radiate 
where they are injected, whereas the lower-energy, radio- 
emitting particles suffer mainly adiabatic losses and radiate 
over the outer parts of the jet, as we have calculated. A 
quantitative study of the acceleration and energy-loss pro- 
cesses over the entire energy spectrum is outside the scope 
of the present paper. 



6 SUMMARY AND CONCLUSIONS 

We have investigated the hypothesis that the radio jets in 
3C 31 can be modelled as adiabatic, decelerating, relativis- 
ts flows. Our technique has several advantages over previous 
work in: modelling linear polarization as well as total inten- 
sity; including the effects of velocity shear; taking account 
of anisotropic emission in the rest frame and fitting to two- 
dimensional images rather than to longitudinal profiles. 

We have found that optimized adiabatic models give a 
fair description of the observed brightness and polarization 
distributions in the outer parts of the modelled region. Their 
fit to the data is inferior to that of the free models of LB, but 
is obtained with many fewer free parameters. The adiabatic 
models cannot describe the inner or flaring regions of the 
jets in 3C31, however: the predicted distributions of total 
intensity and linear polarization are inconsistent with those 
observed. In the innermost region, the jets are clearly non- 
adiabatic unless the emission comes primarily from a part of 
the jet volume which is not expanding with distance from the 
nucleus (we do not resolve the region transversely). In the 
flaring region, a much higher initial velocity is required than 
is allowed by our measurements of jet/counter-jet sidedness 
ratio. The emission in this region is well resolved, and comes 
from the entire jet volume. 

We have shown that a modified adiabatic model can 
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Figure 11. Profiles of X-ray and radio emission and pressure 
along the axis of the main jet. (a) Flux-densit y profiles (on ar- 
bitrary scales) derived as in Fig. 3 of Hardcastl e et alj i2002t) by 
integration in boxes of size 0.3 X 2 arcsec 2 (short axis along the 
jet) on images of resolution 0.6 arcsec FWHM. Points: observed 
X-ray emission (Hardcastle et al. 2002J); dashed line: observed ra- 
dio emission (LB); full line: predicted X-ray emission profile for 
the flaring region alone, as described in the text. Note that the 
radio profile is derived from an image at a lower resolution than 
that used to make the equivalent plot in lHardcastle et al.l (2002) 
in order to match the Chandra point-spread function as closely as 
possible, (b) Profiles of the internal (full), external (dashed) and 
synchrotron minimum (dotted) pres sures from our conservation- 
law analysis lLaing & Bridle 2 0021J . Fig. 4, but projected onto 
the plane of the sky). The vertical dotted lines mark the extent 
of the flaring region as defined by LB and in Section l3.2l 
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still be fitted to the total intensity in this region if we 
add distributed injection of relativistic particles which then 
evolve adiabatically; the region where these particles must 
be injected is also one where there is independent evidence 
for recent particle accelera tion from the detectio n of X- 
ray synchrotron radiation jHardcastle et alJ 12002^ and of 
a local over-pressure in the jet from dynamical arguments 
jLaing fc Bridlell2002bh . 

While the total intensity distributions in the jet and 
counter-jet of 3C 31 can be reproduced satisfactorily by de- 
celerating adiabatic jet models with particle injection in the 
flaring region, the polarization data (degree of linear polar- 
ization and apparent magnetic field direction) cannot. The 
apparent magnetic field configuration contains several fea- 
tures that are qualitatively incompatible with adiabatic evo- 
lution of all the field components in laminar-flow models, 
even in the presence of a velocity shear. We infer that the 
departures from adiabatic conditions in 3C 31 must also in- 
clude deviations either from the flux-freezing description of 
the magnetic fields or from a laminar, axisymmetric velocity 
field. The free models of LB achieved better consistency with 
the observed polarization structure by allowing the field to 
become roughly isotropic at the edges of the jet in the flar- 
ing region. This could be achieved by a turbulent velocity 
component even if the field is frozen into the flow, in which 
case the resulting shear would also contribute to the en- 
hancement in emissivity required to fit the adiabatic models 
in the flaring region. Processes such as dynamo action or 
field-line reconnection might also be significant. 

More detailed analysis of the particle injection process 
inferred for the flaring region would benefit from more sensi- 
tive X-ray, optical and infra-red data and higher-resolution 
radio observations of this region, i.e. from deeper exposures 
with Chandra and HST and from deep higher-resolution po- 
larimetry with the EVLA. Improved knowledge of the in- 
tensity and apparent magnetic field structures within this 
region might also assist development of models for the mag- 
netic field microphysics in this region of the 3C 31 jet, which 
all of our analysis suggests is crucial in determining the de- 
celeration dynamics of the jet, and its brightness and polar- 
ization properties further from the nucleus. 
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